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The null-timelike initial-boundary value problem for a hyperbolic system of equations consists 
of the evolution of data given on an initial characteristic surface and on a timelike worldtube to 
produce a solution in the exterior of the worldtube. We establish the well-posedness of this problem 
for the evolution of a quasilinear scalar wave by means of energy estimates. The treatment is given in 
characteristic coordinates and thus provides a guide for developing stable finite difference algorithms. 
'^—^ • A new technique underlying the approach has potential application to other characteristic initial- 

j boundary value problems. 
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I. INTRODUCTION 



The use of null hypersurfaces as coordinates to describe gravitational waves, as introduced by Bondi [l], was 
key to the understanding and geometric treatment of gravitational waves in the full nonlinear context of general 
relativity di]. In one version of the associated characteristic initial-boundary value problem for Einstein's equations, 
boundary data is given on a timelike worldtube and on an initial outgoing null hypersurface The physical 
picture underlying this null-timelike problem is that the worldtube data represent the outgoing gravitational radiation 
emanating from interior matter sources, while ingoing radiation incident on the system is represented by the initial 
null data. This problem has been developed into a Cauchy-characteristic matching scheme in which the worldtube 
data is supplied by a Cauchy evolution of the interior sources [f|. See @ for a review. Cauchy-characteristic matching 
has been implemented as a numerical evolution code in which the Bondi news function describing the radiation is 
calculated at future null infinity using a finite numerical grid obtained by Penrose compactification 3j. Although 
characteristic evolution codes have successfully simulated many null-timelike problems [6[ and have recently been 
applied to extract the radiation from the inspiral and merger of a binary black hole Q, the well-posedness of the 
null-timelike problem for the Einstein equations has not yet been established. The characteristic formulation of the 
Einstein equations implies that certain variables associated with the radiation satisfy a wave equation. Consequently, 
a necessary condition for the well-posedness of the gravitational problem is that the corresponding problem for the 
t-H . quasilinear wave equation be well-posed. In this paper, as a first step toward treating the gravitational case, we show 
that the quasilinear null-timelike problem for a scalar wave propagating on a curved space background is well posed. 

The characteristic initial value problem did not receive much attention before its importance in general relativity 
was recognized. Historically, the development of computational physics has focused on hydrodynamics, where the 
characteristics typically do not define useful coordinate surfaces and there is no generic outer boundary behavior 
comparable to null infinity. The simplest problem for which the characteristic approach is useful is the Minkowski 
space wave equation, which is satisfied by the components of the fundamental special relativistic fields. Progress 
on the null-timelike problem traces back to Duff 8], where existence and uniqueness was shown for the linear wave 
equation with analytic coefficients and analytic data. Existence and uniqueness was later extended to the C°° case 
of the linear wave equation on an asymptotically flat curved space background by Friedlander [§, [l(| ■ 

The demonstration of well-posedness of the quasilinear boundary problem, i.e. the continuous dependence of the 
solution on the data, depends upon establishing estimates on the derivatives for the linearized problem. This requires 
considering generic lower differential order terms fTTI j . Well-posedness depends crucially on the stability of the problem 
against such lower order perturbations. Otherwise, one cannot localize the problem and use the principle of frozen 
coefficients. 

Partial results estimating the derivatives for characteristic boundary problems were first obtained by Miiller zum 
Hagen and Seifert fl2j. Later Balean carried out a comprehensive study of the differentiability of solutions of the 
null-timelike problem for the flat space wave equation [l3|, [l4| ■ He was able to establish estimates for the derivatives 
tangential to the outgoing null cones but weaker estimates for the time derivatives transverse to the cones had to be 
obtained from a direct integration of the wave equation. The derivatives tangential to the null cone were controlled by 
the derivatives of the data but control of the transverse time derivative required two derivatives of the data. Balean 
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concentrated on the differentiability order of the solution and did not discuss the implications for well-posedness of 
the quasilinear problem. 

Frittelli [15] made the first explicit consideration of well-posedness of the null-timelike problem for the wave equation. 
She adopted the approach of Duff, in which the characteristic formulation of the wave equation is reduced to a canonical 
first order differential form, in close analogue to the symmetric hyperbolic formulation of the Cauchy problem. The 
energy associated with this first order reduction gives estimates for the derivatives of the field tangential to the null 
hypersurfaces. As in Balean's treatment, weaker estimates for the time derivatives were obtained indirectly so that 
well-posedness is not ensured when lower order differential order terms or source terms are included as required for 
the quasilinear she was careful to point out. 

A difficulty underlying the problem can be illustrated in terms of the 1 (spatial)-dimcnsional wave equation 

(af-9l)$ = 0, (1.1) 
where (t, x) are standard space-time coordinates. The conserved energy 

E® = ± J dx((8 i *) 2 + (1.2) 

leads to the well-posedness of the Cauchy problem. In characteristic coordinates (t = t — x, x = i + x), the wave 
equation transforms into 

&0x* = 0. (1.3) 
The conserved energy on the characteristics t — const., 

E(t) = J dx{d x $f, (1.4) 

no longer controls the derivative <9 t <I>. 

Up to now, the only treatment of well-posedness of the characteristic initial value problem valid for the quasilinear 
wave equation has been the work of Rcndall 16] , who considered the double null problem where data is given on a pair 
of intersecting characteristic hypersurfaces. Rendall did not treat the characteristic problem head-on but reduced it to 
a standard Cauchy problem with data on a spacelike hypersurface passing through the intersection of the characteristic 
hypersurfaces. Well-posedness than follows from the classic result for the Cauchy problem. He extended his treatment 
to establish the well-posedness of the double-null formulation of the Einstein gravitational problem. The double null 
problem treated by Rendall is a limiting case of the null-timelike problem considered in this paper. However, Randall's 
approach cannot be applied to the null-timelike problem. Also, the reduction to a Cauchy problem does not provide 
guidance for the development of a stable finite-difference approximation based upon characteristic coordinates. 

Here we consider the null-timelike problem for the quasi-linear wave equation in second differential form in terms of 
characteristic coordinates.. The usual technique for showing that the initial-boundary value problem for a hyperbolic 
system of partial differential equations is well posed is to split the problem into a Cauchy problem and local halfplane 
problems and show that these individual problems are well posed. This works for hyperbolic systems based upon 
a spacelike foliation, in which case signals propagate with finite velocity. Besides the existence and uniqueness of a 
solution, well-posedness implies that the solution depend continuously on the data with respect to an appropriate 
norm. For (|1.1[) . the solutions to the Cauchy problem with compact initial data on t — are square integrable and 
well-posedness can be established using the L 2 norm (|1.2[) . 

However, In characteristic coordinates the 1-dimensional wave equation (|1.3p admits signals traveling in the +x- 
direction with infinite coordinate velocity. In particular, initial data of compact support $(0, x) = f(x) on the 
characteristic t = admits the solution $ = g(t) + f(x), provided that g(0) = 0. Here g(t) represents the profile of 
a wave which travels from past null infinity (x — > — oo) to future null infinity (x —¥ +00). Thus, without a boundary 
condition at past null infinity, there is no unique solution and the Cauchy problem is ill posed. Even with the boundary 
condition $(t, —00) = 0, a source of compact support S(t,x) added to (|1.3p . i.e. 

d t d x $ = S, (1.5) 

produces waves propagating to x — +00 so that although the solution is unique it is still not square integrable. 
On the other hand, consider the modified problem obtained by setting $ = e ax ^, 

d t (d x + a)V = F, V{0,x) = e- ax f(x), a >0 (1.6) 
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where F — e~ ax S. With the boundary condition ^>{t, — oo) = 0, the solutions to (|1.6[) vanish at x = +00 and are 
square integrable. As a result, the Cauchy problem (ll.6[) is well posed with respect to an L 2 norm. For the simple 
example where F = 0, multiplication of (|1.6|) by (2a 1 i l + + ^dt^) and integration by parts gives 

\d t J dx(^d x ^) 2 + 2aH 2 ^ = |y da; ^2(^)5,*- (d f *) 2 ) < 1 1 d^*) 2 . (1.7) 

The resulting inequality 



9 t E < const.£ (1.8) 



for the energy 



(1.9) 



provides the estimates for d x ^ and 'F which are necessary for well-posedness. Estimates for <9f<F, and other higher 
derivatives, follow from applying this approach to the derivatives of (| 1 . 6|) . The approach can be extended to include 
the source term F and other generic lower differential order terms. This allows well-posedness to be extended to the 
case of variable coefficients and, locally in time, to the quasilinear case. 

The 2(spatial)-dimensional model problems considered in Sec. [TT] illustrate how this approach generalizes to the 
multi-dimensional case. We consider the model problems in the modified form analogous to (|1.6[) . By means of this 
technique, the characteristic initial-boundary value problem can again be treated by first considering Cauchy and 
half-plane problems. The demonstration of well-posedness of these model problems presents the underlying ideas in 
a transparent form. 

Our main technique is the use of energy estimates. Although the model problems are treated in the modified form, 
the results can be translated back to the original problem. For example, the modification in going from (|1.5j) to (|1.6|) 
leads to an effective modification of the standard energy for the problem. Rewritten in terms of the original variable 
$ = e aXx &, (|1.9[) corresponds to the energy 

E = - [ dxe- 2ax ({d x <S>) 2 + a 2 <P 2 ] . (1.10) 



2 . 

Thus while the Cauchy problem for (|1.6p is ill posed with respect to the Li norm it is well posed with respect to the 
exponentially weighted norm (|1.10j) . However, rather than modifying the norm, for technical simplicity we deal with 
the modified variable ^ . 

The general arguments presented for our model problems can be applied to a wide range of quasilinear characteristic 
problems. Our motivation for the work here is the application to the null-timelike problem for the quasilinear wave 
equation for a scalar field $ in an asymptotically flat curved space background with source S, 

g ab VaV b $ = S($,d c $,x c ) ) (1.11) 

where the metric g ab and its associated covariant derivative V a are explicitly prescribed functions of ($,:r c ). 
The corresponding flat space wave equation, 



takes the form 



(~d? + dl + d 2 + a|)$ = s, (1.12) 

-(-2d u d r + 9 r 2 )(r$) + -^1— d e (sm 2 6d e $>) + 1 2 flfe = S (1.13) 
r r z sin t> r 2 sin 6 

in null-spherical coordinates (it, r, 9, 0) consisting of a retarded time u — t — r and standard spherical coordinates 
(r,6,cj)). In these coordinates, the Minkowski metric is 

ds 2 = -du 2 - 2dudr + r 2 (dd 2 + sin 2 6d<j) 2 ). (1.14) 

The null-timelike problem consists of determining $ in the region (r > R,u > 0) given data $(u, R, 6, </>) on the 
timelike worldtube r = R and $(0, r, 9, <fi) on the initial null hypersurface u = 0. 

In an asymptotically flat background, the metric (|1 . 14|) generalizes to the Bondi-Sachs form 

g ab dx a dx b = -{e 2fs W - r- 2 h AB W A W B )du 2 - 2e 2f3 dudr - 2h AB W B dudx A + r 2 h AB dx A dx B , (1.15) 
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where x A are angular coordinates such that (u,x A ) = const, along the outgoing null rays. Here the radial coordinate 
r is a surface area coordinate so that the area of the topological spheres (u, r) — const, is 4-7T as measured by Hab- In 
the curved space version of angular coordinates analogous to (|1.14|1 . det (h A B) = sin 2 8. 

In Sec. lIIIl we treat the null-timelike problem for the quasilinear wave equation with Lorentzian metric (|1.15p . 

\-2d u d r + Wd 2 r )(r<5>) + {d r W)d r <$> - ±D A (W A d r <S>) - -^d r (W A D A <S>) + ^D A (e 20 D A <S>) 

= e 2f} S(<S>,d c <5>,x c ), (1.16) 

with initial data $(0,r, x A ) and boundary data &(u, R, x A ). Here Da is the 2-dimensional covariant derivative with 
respect to Hab and the metric coefficients (W, (3,W A , Hab) depend smoothly upon (<E>,w,r, x A ) and the source S 
depends smoothly upon ($, (9 Q <I>, u, r, x A ). 

An essential part of any initial-boundary value problem is the compatibility between the data at the intersection 
between the initial hypersurface and the boundary, i.e. at (u = 0, r = R) in the above case. This compatibility affects 
the differentiability of the resulting solution. In order to avoid difficult issues of analysis, we only give a rigorous 
treatment for the case of smooth initial and boundary data with compact support bounded away from the intersection, 
in which case the solution is C°° locally in time. See the work of Balean [l3|, [ljj for a discussion of the differentiability 
of the solution in the general case. 

We assume that as r — > oo (the approach to null infinity) that the problem reduces to the flat space problem 
(HEI]), so that the coefficients have the asymptotic behavior W = 1 + 0(l/r), ,9 = + 0(l/r), W A = 0(1) and 
Iiab = Qab + 0(l/r), where qAB is the unit sphere metric. The results of Friedlander [l(| then imply that the scalar 
wave falls off as $ ~ $o(w, x A )/r where <E>o is the asymptotic radiation field. 

In Sec. PH we show that the nullcone-worldtube problem for (|1.16|) is well posed subject to the condition that 
S = 0(r~ 3 ) and a positivity condition that the principal part of the wave operator reduces to an elliptic operator in 
the stationary case. Our results are based upon energy estimates obtained by integration by parts with respect to 
the characteristic coordinates. As a result, the analogous finite difference estimates obtained by summation by parts 
provide guidance for the development of a stable numerical evolution algorithm for f| 1 . . 



II. WELL-POSEDNESS OF MODEL CHARACTERISTIC PROBLEMS 



We consider here several model 2(spatial)-dimensional problems which reveal the essential features underlying a 
well posed characteristic initial-boundary value problem. For simplicity of notation, we indicate partial derivatives 
by subscripts, e.g. $t(t, x, y) — <9 t $(£, x, y). Also, we denote the L 2 scalar product and norm over the x, y domain by 
($i,$ 2 ) and ||$|| 2 = ($,$). 

We consider model linear problems with constant coefficients but show that the problems are stable against lower or- 
der perturbations. We also obtain estimates for arbitrarily high derivatives. Thus we can use standard techniques (ll| 
to establish the well-posedness of the corresponding problem with smooth variable coefficients. For the extension to 
the quasilinear case, we require that the coefficients depend smoothly upon the field $ with nonsingular behavior in 
the neighborhood of the initial data. Then well-posedness, locally in time, of the quasilinear problem also follows 
from standard techniques pd| . (See the Appendix of [l7| for details concerning how these standard techniques apply 
to hyperbolic systems in second differential order form.) 

Our goal is to show the well-posedness of the strip problem 

2d> te = ((1 - x)H x ) x + % y + b [{{I - xf$ x ) y + ((1 - x)H y ) x ^j 

in the domain 

< x < 1, -oo < y < oo, t > 

with initial and boundary conditions 

$(0, x, y) = f(x, y), $(t, 0, y) = q(t, y), 

respectively. The method used to show that this problem is well posed applies to the compactified version of the 
null-timelike boundary problem for the wave equation (|1.16[) treated in Sec. IIIII As explained in the Introduction, we 
treat the problem in the modified form obtained by the change of variable $ = e ax ^ , a > 0. 
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A. The Cauchy problem 



We first consider the Cauchy problem 



+ a*) t = % y -2b%, x = (x,y)eR\ t > 0, (2.1) 
y(0,x,y) = f(x,y), (2.2) 

where x and t are both characteristic coordinates. Here a, b are real constants and f(x,y) £ (a smooth function 
with compact support). As explained in the Introduction, we investigate the behavior of square integrable solutions, 
so that \I>(t, ±oo,y) = 0. 

1. The Fourier method 

We first solve the problem by Fourier transform. Let 

f(uj) = —[ e - lQ - s f(x)dxdy, uj = (wi,w 2 )real, 
27T Jb? 

denote the Fourier transform of / and ^(t, lj) the Fourier transform of ^ Then fy(t,U)) is the solution of 

(iwi+a)* t = - (wf + 2biuj 2 ) (2.3) 
*(0,u))=/>), (2.4) 

i.e. 

*t = s*, 

where 

Wo + 26iw9 (wo + 2bi(jj 2 )(a ~ iu>i) . ^ 

s = ~ = ^ ^ . (2.5) 

iwi + a a z + cj^ 

Therefore 



fuJ^Wi — 2a6cji) 



(2.6) 



We now discuss the dependence of the solutions on a, b in detail. 

1) 6 = 0, a > 0. By (|2~6]l . 

w 2 + a 2 

There are no exponentially growing solutions. 

2) b = 0, a = 0. By (|2~^j) . 

5Rs = 0, |9s|^oo for |wi| -> 0, w 2 7^ 0. 
Therefore the solution of (|2.2[) loses all smoothness in time if /(0, 0J2) 7^ 0. 

3) b = 0, a < 0. By (|2~B . 

5Rs +00 for u)\ -> 0, ^2 > 0- 

Thus the problem is ill posed. 
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4) b ^ 0, a > 0. By ([23]) , 

w 2 + a 2 a(w 2 + a 2 ) a 
There is exponential growth but the growth is bounded independently of Co. 

5) b ^ 0, a = 0. By ([^B]) . 

Jis = -. 

Thus there is unbounded exponential growth as uj\ — > 0. The same is true if a < 0. 
We now express our results in a more general setting. 

Definition 2.1. We call the Cauchy problem well posed if, for every f G C^°, i/iere is a unique, smooth, square 
integrable solution and if there is a constant a which does not depend on Co such that 

5Rs < a. 

The problem is ill posed if there is no upper bound a, i.e., there is a sequence QV> such that 

lim 3isj = oo. 

Theorem 2.1. The Cauchy problem i2.2\) is well posed if a > 0. But it is ill posed if a < or a = 0, 6^0. 

2. The energy method 

For the generalization to variable coefficients it is necessary to show that the differential equation (12.2[) is stable 
against lower order perturbations. For this purpose we first apply the energy method to the doubly-characteristic 
Cauchy problem 

(* x + a*) t =^ yy -2b^y-c9 x + d^ t + e^ + F(t,x,y), (2.7) 
<J/(0, x, y) = f(x, y) , -co < x, y < oo , t > 0. 

Here a > 0, b, c, d, e are real constants and F is a forcing (source) term of compact spatial support. 

The term d^t can be absorbed into the left hand side and we obtain (^> x + (a ~ d)^) t . Therefore we neglect this 
term and assume that a is sufficiently large so that a — d > 0. We also neglect the term e^/ because it has no influence 
on the required energy estimates.. Therefore we consider the corresponding Cauchy problem for 

(* x + a*) f = ty yy — 2b^ y — c& x + F. (2.8) 

We now derive an energy estimate. By (j2,8[) . 

(*, V xt ) + o(*, * t ) = -(*«*«) + ^ t ||*|| 2 - % y ) - (*, 26* y + c* K ) + (¥, F). 



Since 



integration by parts gives 



Next, 



^||*|| 2 + ||^ y || 2 = (v^, vj, t ) + < ^II^H 2 + ^||* t || 2 + i(||(v]df + ||F|| 2 ). (2.9) 



(* t , *, t ) + a||* 4 || 2 = -^tll*!/!! 2 ~ 26(* t , * v ) ~ **) + i^t,F). 
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Since 



,•(*,.*,) = < £||* t f + —||^| 



(**,F)< ^||* t || 2 + W, 



we obtain 



Next, 



K„ 1 

T ll* t || 2 + ^II^H 2 < const. (H^ll 2 + \\^ y f + \\Ff) . (2.10) 



Since 0& x ,* w ) = = °, we obtain 

^||*»|| a < const. (H^ll 2 + ||*J 2 + ||F|| 2 ) + |||* t || 2 . (2.11) 
Adding (|2.9p — (|2.1ip gives the energy estimate 

y ||tt t || 2 + ||^f + l dt (ll^H 2 + H^ll 2 + a \\n 2 ) < const. (II^H 2 + II^H 2 + ||vp|| 2 + ||F|| 2 ) . (2.12) 
We have proved: 

Theorem 2.2. The Cauchy problem {2.7}) is well posed with respect to the L2 norm if (a — d) > 0. There is an 
energy estimate. Also, the problem is stable against lower order perturbations. In addition, estimates for the higher 
derivatives of $ follow from the equations obtained by differentiating pT 



We now consider the Cauchy problem 

(* B + o*) t = + * w + F(t, x, y), (2.13) 
®(0,x,y) = f{x,y), ~oo<x,y<oo, i > 0, 

where a; is a characteristic coordinate but t is timelike. We again derive an energy estimate. 
We have 

(*,* ri + a* t ) - -(^,* t ) + |9 t ||*|| 2 

= -(i|vM 2 + nMy 2 ) + oi/,n 



-d t \\n 2 + n^ir + n*vir = *o + (*,*■) < dwr + + (%n (2.14) 



Next. 



i.e. 



(*t, + o* t ) = (* tl *,rt) + a||* 4 || 2 = a||* t || 2 

= -±0 t (||**|| 2 + ||ttj 2 ) + (* t ,n 

«ll* t || 2 + ^(H^H 2 + ll*,H 2 ) = (*t,F) < ^ t \\ 2 + -\\F\\ 2 . (2.15) 

Combining (|2.14|) and (|2.15p as before, we obtain the desired estimate 

^||vl/,|| 2 + ||^|| 2 + ||*„f + ldt (H^ll 2 + H^ll 2 + a||*|| 2 ) < const. (II^H 2 + ||^|| 2 + ||*|| 2 + ||F|| 2 ) . (2.16) 

Remark . As before, we can add a general lower order expression and still obtain the estimate. Also, we can estimate 
all derivatives. 
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B. The half-plane problem 

We now apply the energy method to the double-null halfplane problem for 

(#3 + a*)t = ^ yy + 2b^ y - cl> x + F, 0<x<oo, -oo<y<oo, t>0, (2.17) 
with initial and boundary data 

*(0,x,y,) = f(x,y), *(*,0,!/) = (2.18) 

and source F(t,x,y) of compact support. 

There are no difficulties to derive the basic estimate (|2.12l) because for the estimates (|2.9l) - (|2.1ip we require only 
that ^(t, 0, y) — 0. To obtain estimates for higher derivatives we have to proceed in the following way. 

We differentiate (|2.17[1 with respect to y. Since ^f y (t,0,y) = 0, we obtain the same problem for ty y and therefore 
we obtain estimates for 

ll* ra l| 2 , ll^ll 2 - 

If we differentiate (12.171) two times with respect to y, we obtain estimates for the third derivatives. The corresponding 
results hold for t-derivatives, e.g. 

\\* t f, \\% t f, ||* xt || 2 . 
Now we differentiate (|2.17[) with respect to x. 

(*ra + a$x)t = ^yyx + R- (2.19) 

Here R consists of source terms and terms which we have already estimated. (|2.19[) gives us 

{^xx, ^ xxt) + a(^ xx , $> xt ) = (V X x, ^ yyx) + {^xx, R). (2.20) 

We obtain 

\d t \\^>xx\\ 2 <\UX + « 2 )II^H 2 + ll^'tll 2 + ll* W ,H 2 + PIP) , 

where we already have estimates for H^tH 2 and H^yy^H 2 - The process can be continued. 

Remark . Inhomogeneous boundary data ^{t, 0,y) = q(t, y) may be treated in the same way through the transfor- 
mation — > ^ — qe~ x and absorbing the boundary data in the source term F. We can also treat the timelike-null 
halfplane problem for i2.13\) in the same way. 

C. The strip problem 

As a prototype of the compactified wave equation considered in Sec. IIII1 we consider the strip problem 

2(* a + o*) t - ((1 - x?^x) x + * yy + ^((1 - x?^x) y + ((1 - x) 2 9 y )^ + F(t, x, y) (2.21) 

for 

< x < 1, -co < y < oo, t > 

with initial and boundary conditions 

*(0, x, y) = f(x, y), *(t, 0, y) = q(t, y). 

Here a > and b, with |6| < 1, are real constants and F is a smooth function. The outer boundary T\ at x = 1 is an 
ingoing characteristic so that no boundary condition is allowed. 
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Since the boundary data at To can be absorbed into the source F, we treat the case q — (see the remark in 
Sec. Ill Bj ). We denote the L 2 norm over I\ by 

||*|| 2 ri = J dy* 2 (t,l,y) 
and the Li norm over the boundary ro at x — by 

||*|| 2 r = J dy* 2 (t,0,y). 

We want to show that there is an energy estimate and that the problem is stable against lower order perturbations. 
We derive the necessary estimates. First, 

2(*, + 2a(¥, * t ) = -2(* X1 * t ) + 9 t ||*|| 2 ri + aa t ||*|| 2 

= - ((1 - (1 ~ ~ W^yf - 26^(1 - x)* x , (1 - auJ^J + F), 

a t ||*|| a r 1 +aflt||*H a + ll(l-a:)* IB || a + ll*i,ll 2 
+26^(1 -a;)*,,., (l-x)* w ^ = 2(# x ,* t ) + (*,F). (2.22) 



Next, 



Next, 



Now, 



i.e. 



Also, 



2(* tl + 2a||* t || 2 = ||* t || 2 ri + 2a||tf t || 2 
±9 t (j|(l - x)* x || 2 + ||^|| 2 + 26((1 - a:)*,, (1 - i)^ + (* t) F). 



(2.23) 



2(* X) **t) + 2a(* X) * t ) = aj^H 2 + 2a(# x , * t ) = (j& x , ((1 - x) 2 * x ) x J + (<J X , 

((1 - ar) 2 * x )^ + fcf * x , ((1 - x) 2 Vy) x ) + (* B) F). (2.24) 



# x , ((1 - x) 2 f x ) x ) = -(^ x ,2(l-x)* x J + (* x , (1 - a:) 2 *** 

= -(tf^l-a:)**)- f((l-x) 2 * x ) X) * x ) -||* x || 2 ro 



*,.((] -./■)-*, !, ) = -(*,.(! -.rif, ) -i||* a || 2 r - 



^..((l-*) 2 *,)^ = -fc(«l-£) 2 tf x ) B ,tf x ) =0, 

6^ x ,((l-z) 2 <ig x J = -2&f* x ,(l-a;)tt s J + b(* x ,(l - x) 2 V xy 
= -2b( (1 - x)**,* X 
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Therefore (|2.24l) becomes 

d t \\<f x \\ 2 + (^(l-x)*^ +^ x \\ 2 ra + I||^ !/ || 2 ri +2b((l-x)* x ,* v S J = -2a(y x ,V t ) + {t> X: F). (2.25) 

All the boundary terms have the right sign to enhance the estimates. Therefore we ignore them. Adding the 
simplified estimates (|232l) . ([2^23]) . (|2~25j) gives 



d t {a\\n 2 + ||#x|| 2 + \Qj + Q+ (* X) (1 - x)^ a 
-2h[ (1 - x)y x ,Vy^ + 2(1 - a)(# x , * t ) - 2a||* t || 2 + (* + * t + F) 

<constY||vM 2 + ll*J 2 + ll*l| 2 + ll^l| 2 \ (2 .26) 



where 

Q = ||(1- x)V x \\ 2 + ||*J 2 + 26^(1 - i)* x , (1 - x)*. 
Since |6| < 1, there is a J > such that 

Q> 5(||(i -z)vM 2 + ||vM 2 ). 

Therefore (|2.26l) gives us an energy estimate. We shall now prove that the problem is stable against lower order 
perturbations. We add an expression 

P = A$ x + B% + C** t + 
to (j2~2Ti Then the estimates for (|2~2"2"j) . ((2~2Tt)) and (EHM)) will be changed by lower order terms 

AV X ) + (*, + C^ t ) + £>*) 

A* x ) + (¥ x , BV y ) + (¥ a) C¥ t ) + 

Clearly, there is an energy estimate, provided we choose 2a > \C\, Thus the strip problem (|2.21[) is well posed. 
Now we start with 

2$ xt = ((1 - x)H x ) x + % y + b{(l- x)H x ) y + b ((1 - x)H y ) x + S(t, x, y). (2.27) 
We make the change of variables 



$ x = e ax V x + ae ax *, $ xx = e ox * xx + 2ae ox * x + a 2 e^% 
and set F = e~ ax S. Then we obtain (|2.21[) which is modified by 1Z 

2(* x + o*) t = ((1 - x) 2 *^ + ^ + b (((1 - x) 2 * x ) y + ((1 - x)H y ) x ) +F + TZ. (2.28) 
Here 7Z consists of lower order terms, 

K = 2a(l - x) 2 f x + (a 2 (l - x) 2 - 2a(l - x)) * + 2o6(l - x) 2 ^. 

Since (|2.21l) is stable against lower order terms there is an energy estimate for (|2.28[) . In the same way as in Sec. Ill Bl 
we can estimate the higher derivatives. This allows us to extend well-posedness to the variable coefficient problem 
and, locally in time, to the quasilinear problem. 
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III. THE QUASILINEAR WAVE EQUATION ON AN ASYMPTOTICALLY FLAT BACKGROUND 

We now treat the null-timelike initial-boundary problem (|1.16[) for the quasilinear wave equation. We compactify 
the domain R < r < oo by the transformation x = 1 — R/r to obtain a strip problem < x < 1 with future null 
infinity X + at the boundary x = 1. In terms of the rescaled field $ = r$, the wave equation transforms into 

2d u d x $ - R- x d x (W(l - x) 2 d x $) + (1 - x)R- 1 {d x W)^ 
+R- 2 D A ((1 ~ x) 2 W A d x §) + R- 2 d x ((l - x) 2 W A D A §) - R~ 2 (l - x)(D A W A )$ 

-R- 1 D A {e 2f} D A ^>) = ~r 3 i?~ 1 e 2 ' 9 5. (3.1) 

Here we use the 2-metric 1i A b and its inverse h AB to raise and lower indices of tensor fields on the spacclike (u — 
const., r = const.) spherical cross-sections . Up to lower order terms, (|3.1|) is a 3(spatial)-dimensional version of (|2.27j) 
where the y-coordinate has been replaced by the x^-coordinate on the spherical cross-sections and the t-coordinatc 
has been replaced by the u-coordinate. In order for our treatment to apply to the quasilinear case, we assume that 
the metric coefficients (W, f3, W A , h A s) depend smoothly upon ($, u, r, x A ) and that the source S depends smoothly 
upon (<f>, <9 a $, it, r, x A ), with non-singular Lorentzian geometry in the neighborhood of the initial data. 

We treat the modified problem resulting from the transformation <t = e ax ^. The same argument used in Sec. Ill Cl 
shows that this problem is stable with respect to lower order terms. We ignore these terms and thus obtain the strip 
problem 

2d v (d x y + a*) = R~ l d x {W{l - x) 2 d x t>) - R- 2 D A {{1 - x) 2 W A d x t>) - R-' 2 d x {{\ - x) 2 W A D A ^) 

+R- 1 D A (e 20 D A ^) +F, *(0,x,x A ) = /, V(u,0,x A ) = q , (3.2) 

where F = —r 3 R~ 1 e 2/3 e ax S. In order to treat (|3.2p . we require that the physical space source has asymptotic behavior 
S = 0(r~ 3 ) so that F is square integrable over the strip. No boundary condition is allowed at the outer boundary 
Tx at x = 1 since T + is an ingoing characteristic surface. 

We obtain the required estimates for (|3.2p by the same method used in Sec. Ill CI The data at the inner boundary 
Fo at x = may be absorbed into F so it suffices to treat the case q = 0. We define the inner product 

(*i,* 2 ) = J dx <j> dw*i* 2 

and L2 norm H^H 2 = ("J, where dw is the area element on the unit sphere. We write 

\\V A \\ 2 = (V A ,V A ) = (h AB V A ,V B ). 

Since the spherical cross-sections are spacelike, their intrinsic 2-metric h A B is positive definite so that || V A \\ serves as 
an Li norm for the angular components. We also need a metric norm for spacelike 3- vectors. In the standard Cauchy 
problem this is supplied by the intrinsic 3-metric of the spacelike Cauchy hypersurfaces. Since the characteristic 
hypersurfaces have a degenerate 3-metric, we take a different approach. We use the projection operator ir£ = S£—t a dbU, 
where t a d a = d u , to define a 3-metric r y ab = ■K < ^TT b d g cd . For the Bondi-Sachs metric (|1 . 15|) . the resulting components in 
the (u,r, x A ) coordinates are 

y™ = 0, 7 rr = e- 2 PW, Y A = -e- 2 ?r- 2 W A , 7 AB = r - 2 h AB . 

Denoting x % — (r, X A ), this implies that 7 l - J — > e 4 - 7 as r — > oo, where e 4 - 7 is the Euclidean 3-metric expressed in standard 
spherical coordinates. In the compactified coordinates, x l = (x, x A ) it is more useful to deal with the rescaled 3-metric 
j ab = e 2 P r 2 j ab which has components 

f™ = 0, r x = (1 - x) 2 W, y cA = -R- 1 {l-xfW A , ^ AB = e 2(s h AB . 

We then define 

mf = (rv i ,v j ) 

which serves as an norm for the x % components. Thus 

\\&m 2 = ||W X / 2 (1 - x)d x V\\ 2 + \\ePD A *f - 2R- 1 ((1 - x)d x V, (1 - x)W A D A *) . (3.3) 
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We also define the corresponding inner products and norms on the boundaries, e.g. 

(*i,* 2 )r = j> dw*i* 2 , ||*||r = (*,*)r- 

Because the original radial coordinate r was a surface area coordinate, the 2-metric Hab of the spherical cross- 
sections has determinant det{h A B) = det{q A s), where q A B is the unit sphere metric. Consequently, 

D A D A ^ 2 ) = -(D A y u D A y 2 ) (3.4) 

and 

(V A ,D A 9) = -{D A V A ,9) (3.5) 

where V A {u 1 x 1 x A ) is any smooth vector field on the spherical cross-sections. These identities allow the necessary 
integration by parts. 

We derive the required estimates by freezing the dependence of the metric coefficients on (^,u, x) but we retain 
their dependence on x A so that W A and h A B remain smooth vector and tensor fields on the spherical cross-sections. 
We follow the procedure in Sec. Ill CI First, 

2(%a u a a 9) + 2a(*,a u *) = -2(d x *,a u v) + 9 u ||*|| 2 ri + ad u \n 2 

= -R- 1 \W(1 - x)d x V, (1 - x)d x *j - R-^Da^W 2 + 2Rr 2 ^(1 - x)W A D A V, (1 - x)d x vj + (^,F), 

duWn 2 ^ + ad u \\n 2 + R-'Wd^f = 2(9,*, du V) + (tf ,F). (3.6) 

Next, 

2(d u y,d u d x *) + 2a\\d u n 2 = R*|| 2 ri +2a||d„*|| 2 
= -^R-'du (\\W^ 2 (1 - x)d x n 2 + \\e p D A ^f - 2R-\{1 - x)d x V, (1 - x)W A D A ^) + {d u V,F) , 



so that 



Next, 



|3 u *|| 2 ri +2a\\d u n 2 + \R- 1 du\\d i n 2 = (du'S>,F). (3.7) 



2(d x V, d u d x V) + 2a(d x V, d u V) = d u \\d x n 2 + 2a(d x V, d u ^) 
= R- 1 (d x 9, d x (W(l - x) 2 d x *)) + R- 1 D A {e^D A *) 



-R- 2 [y x ,D A {{l-xfW A d x V)) -R- 2 (d x *,d x ((l-x) 2 W A D A y) )+(d s %F). 
As shown in in Sec. Ill CI 



Also, 



d x *, d x (W(l -x) 2 d x *) \ =-(^%W(l-x)d x ^J-^\\W 1 ^d x 9\\ 2 r . 



(d x % D A (e 2p D A ^)) = -I|| e ^ A *|| 2 ri (3.9) 



and 

(d„*,D A ((l-x) 2 W A d s *)) + [d a %d x ((l-x) 2 W A D A V) I =-2 1 (l-x)d m *,W A D A 9 
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Therefore (13.81) becomes 

0„||0x*|| a + R- 1 (d x V, W(l - x)d x *j + iR-^W^d^fr, + ii?- 1 ||e' 3 J D yl *|| 2 r 1 

-2iT 2 ( (1 - x)d x *,W A D A y) =-2a(d x V,d u *) + (d x y,F). (3.10) 



As before, the boundary terms have the right sign to enhance the estimates so that we can ignore them. Adding 
the simplified estimates (|3.6p . (I3.7[) . (|3.10j) gives 

d u (a\\n 2 + \\d x n 2 + l^W^n 2 ^) + ir 1 !^ 2 + r- 1 (a**, w(i - x)a x *J 

= 2R- 2 ^(1 - x)d a V, W A D A *J + 2(1 - o)^*, flu*) - 2<z||d„*|| 2 + (* + d u ^> + d x *, F) 

< const/ 1|^|| 2 + H^a^II 2 + ||*|| 2 + \\F\A • (3.11) 



Therefore (|3.11[) gives us an energy estimate provided that the 3-metric 7 y has (+ + +) signature, so that is 
a norm for the gradient d^f — (d r ^, 8a^)- This is equivalent to the requirement that the principle part of the wave 
operator reduce to an elliptic operator in the stationary case where the it-derivatives vanish. Since 7 IJ is asymptotic 
to the Euclidean metric as r — > oo, this positive-definite condition is satisfied throughout some exterior domain. 

Estimates for the higher derivatives of ^ and stability against lower order perturbations follow from the same 
arguments given in Sec. [TTJ This establishes the well-posedness of the worldtube-nullcone problem for the case of 
smooth variable coefficients. The extension of well-posedness, locally in time, for the quasilinear case then follows 
from the standard techniques referred to in Sec. [TT] . 

For a mass M Schwarzschild geometry, j rr — e~ 2/3 W = 1 — 2M/r so that positive-definiteness of the 3-metric ^ 
breaks down at r = 2M where the worldtube becomes null. In this limiting case of the double-null problem, the 
<9«||9 X 4'|| 2 term in (I3.11[) suffices to provide the required estimate. However, for R < 2M the "worldtube" is spacelike 
and the problem must be treated differently. 
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The Well-posedness of the Null-Timelike Boundary Problem for Quasilinear Waves 
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The null-timelike initial-boundary value problem for a hyperbolic system of equations consists 
of the evolution of data given on an initial characteristic surface and on a timelike worldtube to 
produce a solution in the exterior of the worldtube. We establish the well-posedness of this problem 
for the evolution of a quasilinear scalar wave by means of energy estimates. The treatment is given in 
characteristic coordinates and thus provides a guide for developing stable finite difference algorithms. 
'^—^ • A new technique underlying the approach has potential application to other characteristic initial- 

j boundary value problems. 
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I. INTRODUCTION 



The use of null hypersurfaces as coordinates to describe gravitational waves, as introduced by Bondi [l], was 
key to the understanding and geometric treatment of gravitational waves in the full nonlinear context of general 
relativity di]. In one version of the associated characteristic initial-boundary value problem for Einstein's equations, 
boundary data is given on a timelike worldtube and on an initial outgoing null hypersurface The physical 
picture underlying this null-timelike problem is that the worldtube data represent the outgoing gravitational radiation 
emanating from interior matter sources, while ingoing radiation incident on the system is represented by the initial 
null data. This problem has been developed into a Cauchy-characteristic matching scheme in which the worldtube 
data is supplied by a Cauchy evolution of the interior sources [f|. See @ for a review. Cauchy-characteristic matching 
has been implemented as a numerical evolution code in which the Bondi news function describing the radiation is 
calculated at future null infinity using a finite numerical grid obtained by Penrose compactification 3j. Although 
characteristic evolution codes have successfully simulated many null-timelike problems [6[ and have recently been 
applied to extract the radiation from the inspiral and merger of a binary black hole Q, the well-posedness of the 
null-timelike problem for the Einstein equations has not yet been established. The characteristic formulation of the 
Einstein equations implies that certain variables associated with the radiation satisfy a wave equation. Consequently, 
a necessary condition for the well-posedness of the gravitational problem is that the corresponding problem for the 
t-H . quasilinear wave equation be well-posed. In this paper, as a first step toward treating the gravitational case, we show 
that the quasilinear null-timelike problem for a scalar wave propagating on a curved space background is well posed. 

The characteristic initial value problem did not receive much attention before its importance in general relativity 
was recognized. Historically, the development of computational physics has focused on hydrodynamics, where the 
characteristics typically do not define useful coordinate surfaces and there is no generic outer boundary behavior 
comparable to null infinity. The simplest problem for which the characteristic approach is useful is the Minkowski 
space wave equation, which is satisfied by the components of the fundamental special relativistic fields. Progress 
on the null-timelike problem traces back to Duff 8], where existence and uniqueness was shown for the linear wave 
equation with analytic coefficients and analytic data. Existence and uniqueness was later extended to the C°° case 
of the linear wave equation on an asymptotically flat curved space background by Friedlander [§, [l(| ■ 

The demonstration of well-posedness of the quasilinear boundary problem, i.e. the continuous dependence of the 
solution on the data, depends upon establishing estimates on the derivatives for the linearized problem. This requires 
considering generic lower differential order terms fTTI j . Well-posedness depends crucially on the stability of the problem 
against such lower order perturbations. Otherwise, one cannot localize the problem and use the principle of frozen 
coefficients. 

Partial results estimating the derivatives for characteristic boundary problems were first obtained by Miiller zum 
Hagen and Seifert fl2j. Later Balean carried out a comprehensive study of the differentiability of solutions of the 
null-timelike problem for the fiat space wave equation [l3|, [l4| ■ He was able to establish estimates for the derivatives 
tangential to the outgoing null cones but weaker estimates for the time derivatives transverse to the cones had to be 
obtained from a direct integration of the wave equation. The derivatives tangential to the null cone were controlled by 
the derivatives of the data but control of the transverse time derivative required two derivatives of the data. Balean 
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concentrated on the differentiability order of the solution and did not discuss the implications for well-posedness of 
the quasilinear problem. 

Frittelli [15] made the first explicit consideration of well-posedness of the null-timelike problem for the wave equation. 
She adopted the approach of Duff, in which the characteristic formulation of the wave equation is reduced to a canonical 
first order differential form, in close analogue to the symmetric hyperbolic formulation of the Cauchy problem. The 
energy associated with this first order reduction gives estimates for the derivatives of the field tangential to the null 
hypersurfaces. As in Balean's treatment, weaker estimates for the time derivatives were obtained indirectly so that 
well-posedness is not ensured when lower order differential order terms or source terms are included as required for 
the quasilinear she was careful to point out. 

A difficulty underlying the problem can be illustrated in terms of the 1 (spatial)-dimcnsional wave equation 

(af-9l)$ = 0, (1.1) 
where (t, x) are standard space-time coordinates. The conserved energy 

E® = ± J dx((8 i *) 2 + (1.2) 

leads to the well-posedness of the Cauchy problem. In characteristic coordinates (t = t — x, x = i + x), the wave 
equation transforms into 

&0x* = 0. (1.3) 
The conserved energy on the characteristics t — const., 

E(t) = J dx{d x $f, (1.4) 

no longer controls the derivative <9 t <I>. 

The first proof of well-posedness of the characteristic initial value problem valid for the quasilinear wave equation 
has been the work of Rcndall |16l | , who considered the double null problem where data is given on a pair of intersecting 
characteristic hypersurfaces. Rendall did not treat the characteristic problem head-on but reduced it to a standard 
Cauchy problem with data on a spacelike hypersurface passing through the intersection of the characteristic hyper- 
surfaces. Well-posedness than follows from the classic result for the Cauchy problem. He extended his treatment to 
establish the well-posedness of the double-null formulation of the Einstein gravitational problem. The double null 
problem treated by Rendall is a limiting case of the null-timelike problem considered in this paper. However, Rendall's 
approach cannot be applied to the null-timelike problem. Also, the reduction to a Cauchy problem does not provide 
guidance for the development of a stable finite-difference approximation based upon characteristic coordinates. 

Another limiting case of the null-timelike problem is the Cauchy problem on a characteristic cone, corresponding 
to the limit in which the timelike worldtube has shrunk to a nonsingular worldline. This problem is difficult to treat 
in characteristic coordinates because of their singular nature at the vertex of the cone. However, Choquet-Bruhat, 
Chrusciel and Martm-Garcfa have been able to establish the existence of solutions to this problem, for both the scalar 
and gravitational case, by treating it in harmonic coordinates adapted to the null cones |17l . [l8j . 

Here we consider the null-timelike problem for the quasi-linear wave equation in second differential form in terms of 
characteristic coordinates.. The usual technique for showing that the initial-boundary value problem for a hyperbolic 
system of partial differential equations is well posed is to split the problem into a Cauchy problem and local halfplane 
problems and show that these individual problems are well posed. This works for hyperbolic systems based upon 
a spacelike foliation, in which case signals propagate with finite velocity. Besides the existence and uniqueness of a 
solution, well-posedness implies that the solution depend continuously on the data with respect to an appropriate 
norm. For (|1.1[) . the solutions to the Cauchy problem with compact initial data on t — are square integrable and 
well-posedness can be established using the Li norm (|1.2[) . 

However, In characteristic coordinates the 1-dimensional wave equation (|1.3p admits signals traveling in the +x- 
direction with infinite coordinate velocity. In particular, initial data of compact support $(0, x) — f(x) on the 
characteristic t — admits the solution $ = g(t) + f(x), provided that g(0) — 0. Here g{t) represents the profile of 
a wave which travels from past null infinity (x — > — oo) to future null infinity (x — > +oo). Thus, without a boundary 
condition at past null infinity, there is no unique solution and the Cauchy problem is ill posed. Even with the boundary 
condition $(i, — oo) = 0, a source of compact support S(t, x) added to (II. 3[) . i.e. 



d t d x <!> = S, 



(1.5) 
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produces waves propagating to x — +00 so that although the solution is unique it is still not square integrable. 
On the other hand, consider the modified problem obtained by setting <I> = e ax ^, 

d t (d x +a)V = F, V(0,x) = e- ax f{x), a > (1.6) 

where F — e~ ax S. With the boundary condition \I>(t, —00) = 0, the solutions to (|1.6[) vanish at x = +00 and are 
square integrable. As a result, the Cauchy problem (jl.6|) is well posed with respect to an L 2 norm. For the simple 
example where F = 0, multiplication of (11.61) by (2a\& + d x ^ + ^9^) and integration by parts gives 



\d t J dx((d x *f + 2a 2 ^ = ~j dx(2(d t V)d x * - (dt*) 2 ) <~J dx{d x *f 



(1.7) 



The resulting inequality 



d t E < const.E (1.8) 



for the energy 

E =\ J dx(^(d x ^) 2 + 2a 2 ¥ 



(1.9) 



provides the estimates for d x ^> and ^ which are necessary for well-posedness. Estimates for dt*&, and other higher 
derivatives, follow from applying this approach to the derivatives of (jT76j) . The approach can be extended to include 
the source term F and other generic lower differential order terms. This allows well-posedness to be extended to the 
case of variable coefficients and, locally in time, to the quasilinear case. 

The 2(spatial)-dimensional model problems considered in Sec. [TT| illustrate how this approach generalizes to the 
multi-dimensional case. We consider the model problems in the modified form analogous to (|1.6p . By means of this 
technique, the characteristic initial-boundary value problem can again be treated by first considering Cauchy and 
half-plane problems. The demonstration of well-posedness of these model problems presents the underlying ideas in 
a transparent form. 

Our main technique is the use of energy estimates. Although the model problems are treated in the modified form, 
the results can be translated back to the original problem. For example, the modification in going from (|1.5I) to (|1.6[) 
leads to an effective modification of the standard energy for the problem. Rewritten in terms of the original variable 
$ = e ax ^, (|1.9p corresponds to the energy 

E = - [ dxe- 2ax ({d x <5>) 2 + a 2 $A . (1.10) 



2 . 

Thus while the Cauchy problem for (|1.6p is ill posed with respect to the Li norm it is well posed with respect to the 
exponentially weighted norm (jl.lOp . However, rather than modifying the norm, for technical simplicity we deal with 
the modified variable ^ . 

The general arguments presented for our model problems can be applied to a wide range of quasilinear characteristic 
problems. Our motivation for the work here is the application to the null-timelike problem for the quasilinear wave 
equation for a scalar field $ in an asymptotically flat curved space background with source S 1 , 

g ab V a V b $ = S($,d c $,x c ) ) (1.11) 

where the metric g ab and its associated covariant derivative V a are explicitly prescribed functions of ($,x c ). 
The corresponding flat space wave equation, 

(~d 2 + d 2 +d 2 +d 2 )<S> = S, (1.12) 

takes the form 

-{-2d u d r + d 2 )(r<S>) + —}—d e {sm6d e ®) + 1 2 flgg = S (1.13) 
r r z sin r z sin 

in null-spherical coordinates (it, r, 9, (f>) consisting of a retarded time u = t — r and standard spherical coordinates 
(r,8,(j)). In these coordinates, the Minkowski metric is 

ds 2 = -du 2 - 2dudr + r 2 (d8 2 + sin 2 8dcj) 2 ) . (1.14) 
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The null-timelike problem consists of determining $ in the region (r > R,u > 0) given data $(«, i?, £?, 0) on the 
timelike worldtube r = R and $(0, r, 9, <fi) on the initial null hypersurface u = 0. 

In an asymptotically flat background, the metric (|1.14j) generalizes to the Bondi-Sachs form 

g ab dx a dx b = -{e 2p W - r - 2 h AB W A W B )du 2 - 2e 2p dudr - 2h AB W B dudx A + r 2 h AB dx A dx B , (1.15) 

where x A are angular coordinates such that (u,x A ) = const, along the outgoing null rays. Here the radial coordinate 
r is a surface area coordinate so that the area of the topological spheres (u, r) = const, is 4ir as measured by the 
conformal 2-metric h AB . In the curved space version of angular coordinates analogous to (|1.14|) . det(h AB ) = sin 2 9. 

In Sec. IIII1 we treat the null-timelike problem for the quasilinear wave equation (jl.ll[) with asymptotically flat 
Lorentzian metric (|1 . 15|) . 

\-2d u d r + Wd 2 r ){r<S>) + {d r W)d r <S> - ^D A {W A d r <b)-^d r {W A D A §) + ^D A {e 20 D A >S>) 

= e 2/3 S($,d c $,x c ), (1.16) 
$(0,r,a; A ) = f{r,x A ) , $(u, R, x A ) = q{u, x A ) , R < r < oo , u>0. 

Here D A is the 2-dimensional covariant derivative with respect to h AB and the metric coefficients (W, (3,W A , h AB ) 
depend smoothly upon ($,ii,r, x A ) and the source S depends smoothly upon ($, d a <&, u, r, x A ). 

An essential part of any initial-boundary value problem is the compatibility between the data at the intersection 
between the initial hypersurface and the boundary, i.e. at (u = 0, r — R) in the above case. This compatibility affects 
the differentiability of the resulting solution. In order to avoid difficult issues of analysis, we only give a rigorous 
treatment for the case of smooth initial and boundary data with compact support bounded away from the intersection, 
in which case the solution is C°° locally in time. See the work of Balean [l3|, [lj| for a discussion of the differentiability 
of the solution in the general case. 

We assume that as r -> oo (the approach to null infinity) that the problem reduces to the flat space problem 
(I1.13[) . so that the coefficients have the asymptotic behavior W = 1 + 0(l/r), (3 = + 0(l/r), W A = 0(1) and 
h AB = q AB + 0(l/r), where q AB is the unit sphere metric. The results of Friedlander [l(| then imply that the scalar 
wave falls off as $ ~ $o( M : x A )/r where $o is the asymptotic radiation field. 

In Sec. Jm we establish our Main Theorem: 

The nullcone- worldtube problem (| 1 . 1 6|) is well posed for smooth, compatible initial data f(r,x A ) and boundary data 
q(u,x A ) subject to the conditions that / = 0(r _1 ), S = 0(r~ 3 ) and a positivity condition that the principal part of 
the wave operator reduces to an elliptic operator in the stationary case. 

Our treatment is based upon energy estimates obtained by integration by parts with respect to the characteristic 
coordinates. As a result, the analogous finite difference estimates obtained by summation by parts provide guidance 
for the development of a stable numerical evolution algorithm for (II. 16ft . 



II. WELL-POSEDNESS OF MODEL CHARACTERISTIC PROBLEMS 



We consider here several model 2(spatial)-dimensional problems which reveal the essential features underlying a 
well posed characteristic initial-boundary value problem. For simplicity of notation, we indicate partial derivatives 
by subscripts, e.g. $*(£, x, y) = 9($(i, x, y). Also, we denote the Li scalar product and norm over the x, y domain by 
($i,$a) and ||$|| 2 = (*,$). 

We consider model linear problems with constant coefficients but show that the problems are stable against lower or- 
der perturbations. We also obtain estimates for arbitrarily high derivatives. Thus we can use standard techniques [ll| 
to establish the well-posedness of the corresponding problem with smooth variable coefficients. For the extension to 
the quasilinear case, we require that the coefficients depend smoothly upon the field $ with nonsingular behavior in 
the neighborhood of the initial data. Then well-posedness, locally in time, of the quasilinear problem also follows 
from standard techniques [llT | . (See the Appendix of [l9| for details concerning how these standard techniques apply 
to hyperbolic systems in second differential order form.) 

Our goal is to show the well-posedness of the strip problem 

24> te - ((1 - x) 2 <5> x ) x + +&[((!- x) 2 <5> x ) + ((1 - xf^ y ) x ) 
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in the domain 



< x < 1, -oo < y < oo, £ > 

with initial and boundary conditions 

*(0,a!,|/) = /(a!,tf), $(t,0,i/) = g(t,|/), 

respectively. The method used to show that this problem is well posed applies to the compactified version of the 
null-timelike boundary problem for the wave equation (|1.16|) treated in Sec. Mil As explained in the Introduction, we 
treat the problem in the modified form obtained by the change of variable $ = e ax ^f, a > 0. 



A. The Cauchy problem 



We first consider the Cauchy problem 

+ a*) t = V yy - 2b%, x = (x, y) G i? 2 , t > 0, (2.1) 
*(0,x,y) = f(x,y), (2.2) 

where x and t are both characteristic coordinates. Here a, b are real constants and f(x,y) £ C§° (a smooth function 
with compact support). As explained in the Introduction, we investigate the behavior of square integrable solutions, 
so that ±oo, y) = 0. 



1. The Fourier method 

We first solve the problem by Fourier transform. Let 

f(w) = — e- ia - s f(x)dxdy, Q = (cj u oj 2 ) real, 
27T Jb? 

denote the Fourier transform of / and $(t, u>) the Fourier transform of Then \&(t, u>) is the solution of 

+ a)$ t = - (w| + 2biu 2 ) (2.3) 



i.e. 



*(0,w) = /(w), (2.4) 



where 

u>2 + 2biui2 (w| + 2biu>2){a, — icj\) 
iuji + a a 2 + ujf 

Therefore 



We now discuss the dependence of the solutions on a, b in detail. 
1) 6 = 0, a > 0. By ([2~6]h 



5R S = -4^<0. 



There are no exponentially growing solutions. 



(2.5) 



= (^lZ 2 ^il. (2.6) 
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2) b = 0, a = 0. By (|^6|) . 

5Rs = 0, |9s|^oo for -> 0, w 2 7^ 0. 

Therefore the solution of (12. 2[) loses all smoothness in time if /(0, W2) 7^ 0. 

3) 6 = 0, a < 0. By ([2~6]) . 

Sis — > +00 for aj2 ^ co- 
Thus there is unbounded exponential growth and the problem is ill posed. 

4) b 7^ 0, a > 0. By (l2~6l) . 

a ,_.4K> + *f + *4 <g. 

There is exponential growth but the growth is bounded independently of Co. 

5) 6 ^ 0, a = 0. By (l2~B . 



Thus there is unbounded exponential growth as io\ — > 0. The same is true if a < 0. 
We now express our results in a more general setting. 

Definition 2.1. We call the Cauchy problem well posed if, for every f £ C^°, there is a unique, smooth, square 
integrable solution and if there is a constant a which does not depend on Co such that 

5Rs < a. 

The problem is ill posed if there is no upper bound a, i.e., there is a sequence uw) such that 

lim Sftsj = 00. 

Theorem 2.1. The Cauchy problem \2. 2\) is well posed if a > 0. But it is ill posed if a < or a — 0, b 7^ 0. 

2. The energy method 

For the generalization to variable coefficients it is necessary to show that the differential equation (|2.2j) is stable 
against lower order perturbations. For this purpose we first apply the energy method to the doubly-characteristic 
Cauchy problem 

(^ :c +a^) t =^ yy -2b^ y -c^ x + d^t + e^ + F(t,x,y), (2.7) 
<J/(0, x, y) = f(x, y) , —00 < x, y < 00 , t > 0. 

Here a > 0, b, c, d, e are real constants and F is a forcing (source) term of compact spatial support. 

The term d^t can be absorbed into the left hand side and we obtain (^> x + (a — d)^) t . Therefore we neglect this 
term and assume that a is sufficiently large so that a — d > 0. We also neglect the term e^ because it has no influence 
on the required energy estimates.. Therefore we consider the corresponding Cauchy problem for 

+ o#) ( = V yy - 2bV y - c* x + F. (2.8) 

We now derive an energy estimate. By (|2.8[) . 
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Since 



integration by parts gives 



2 iT _ \/a, T , "\ ^ 2 |KTi ||2 i a |liTi l|2 



« 



Next, 



Since 



^ t ||v]>|| 2 + ||*,|| 2 = (* X) %) + (%F) < l\\* x \\ 2 + l\\*t\\ 2 + l(\\m 2 + \\F\\ 2 )- (2-9) 



(* t) * xt ) + a||*t|| 2 = -\dt\\^ y \\ 2 - 2b(V t , * w ) - c(* tl * B ) + F). 



(^^)<^||* t || 2 + -||F|| 2 , 
8 a 



we obtain 



Next, 



2 Vo V " 8 a 



¥ ll* t || 2 + 2^H*»II 2 < const. (||*M 2 + ||^|| 2 + ||F|| 2 ) . (2.10) 



(**, *xt) + a(*», *t) = *«,) ~ 26* y + c$ x ) + F). 
Since (<&x,* w ) = -(^j/,^) = 0, we obtain 

^||*x|| a < const. (H^ll 2 + ||* B f + ||F|| 2 ) + |||* t || 2 . (2.11) 
Adding (I2.9[) - (I2.11[) gives the energy estimate 

y ll* t || 2 + ||*«f + \dt (||*,|| 2 + ||*,|| 2 + a||*|| 2 ) < const. (||^|| 2 + ||*J 2 + ||M/|| 2 + ||F|| 2 ) . (2.12) 
We have proved: 

Theorem 2.2. The Cauchy problem is well posed with respect to the Li norm if (a — d) > 0. There is an 

energy estimate. Also, the problem is stable against lower order perturbations. In addition, estimates for the higher 
derivatives of \ff follow from the equations obtained by differentiating \2. 7] ) . 

We now consider the Cauchy problem 

(* a + o*)t = * ra + + F(t, x, y), (2.13) 
y(0,x,y) = f(x,y), -oo<x,y<oo, t>0, 

where x is a characteristic coordinate but t is timelike. We again derive an energy estimate. 
We have 

= -(||^ + ||^||2 ) + ( ^ n 
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~d 4 ||*|| 2 + II^H 2 + \\* y f = (* X) * t ) + (*,F) < ^||* t || 2 + h* a f + (2.14) 
z o a 



Next, 



i.e. 



= -±W*ll 2 + M 2 ) + (*t,n 

a||^|| 2 + U(||\M 2 + ||^|| 2 ) = (%,F) < |||* t || 2 + -\\FW 2 . (2.15) 

Combining (|2.14[) and f|2.15D as before, we obtain the desired estimate 

^||* t || 2 + ||^|| 2 + ll^f + \d t (H^ll 2 + ||* y || 2 + a\\n 2 ) < const. (||* T || 2 + ||vl/,J| 2 + ||vl/|| 2 + || F || 2 ) . (2.16) 

Remark . As before, we can add a general lower order expression and still obtain the estimate. Also, we can estimate 
all derivatives. 

B. The half-plane problem 



We now apply the energy method to the double-null halfplane problem for 

+ a*) t = tyyy + 2b^ y - c# x + F , 0<x<oo, -oo<y<oo, t>0, (2.17) 
with initial and boundary data 

*(Q,x,y,) = f(x,y), 9{t,0,v) = (2-18) 

and source F(t,x,y) of compact support. 

There are no difficulties to derive the basic estimate (|2.12p because for the estimates (|2.9p - (|2.11l) we require only 
that ^(t, 0, y) = 0. To obtain estimates for higher derivatives we have to proceed in the following way. 

We differentiate (j2 . 1 7[) with respect to y. Since \& y (t, 0, j/) = 0, we obtain the same problem for ty y and therefore 
we obtain estimates for 

\\%yf, ll^ll 2 - 

If we differentiate (|2.17l) two times with respect to y, we obtain estimates for the third derivatives. The corresponding 
results hold for i-derivatives, e.g. 

||vl/ 4 || 2 , ||^|| 2 , ||* xt || 2 . 
Now we differentiate (I2.17[) with respect to x. 

(y xx + a$ x ) t = Vyy X + R. (2.19) 
Here R consists of source terms and terms which we have already estimated. (I2.19[) gives us 

(Vxx, Vxxt) + *xt) = (*«,* W x) + (Vxx, R)- (2-20) 

We obtain 



^||*»*|| 2 < 77 (1 + 2 )||*x*|| a + ll^tll 2 + ll* ra ,|| 2 + Pll 2 



ld t \\* xx f<l 

where we already have estimates for H^xtll 2 and ||^/ ay2: || 2 . The process can be continued. 

Remark . Inhomogeneous boundary data ^(t,0,y) = q{t,y) may be treated in the same way through the transfor- 
mation ^ — > ^ — qe~ x and absorbing the boundary data in the source term F. We can also treat the timelike-null 
halfplane problem for H2.13\) in the same way. 
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C. The strip problem 

As a prototype of the compactified wave equation considered in Sec. IIIIl we consider the strip problem 

2(* a + o*) t = ((1 - x)H x ) x + * y , + b(({l - x)H x ) y + ((1 - x) 2 * y )^j + F(t,x,y) (2.21) 

for 

< x < 1, -co < y < oo, £ > 

with initial and boundary conditions 

*(0, x, y) = /(x, j/), *(t, 0, y) = g(t, y). 

Here a > and 6, with |6| < 1, are real constants and F is a smooth function. The outer boundary ri at x — 1 is an 
ingoing characteristic so that no boundary condition is allowed. 

Since the boundary data at Tq can be absorbed into the source F, we treat the case q = (see the remark in 
Sec. Ill Bj) . We denote the L2 norm over Ti by 

||*|| 2 ri - J dy&(t,l,y) 
and the L 2 norm over the boundary r at x = by 

||*|| 2 r = J dy* 2 (t,0,y). 

We want to show that there is an energy estimate and that the problem is stable against lower order perturbations. 
We derive the necessary estimates. First, 

2(*, + 2o(*, * t ) = -2(V X , * t ) + 9 t ||*|| 2 ri + ad t \\V\\ 2 

= - ((1 - x)9 x , (1 - a:)**) ~ \\*J 2 - 26^(1 - x)*„ (1 - x)^ + F), 

^ll^llV+aftll^f + IKl-^lP + II^H 2 
+2&((1 - z)tf x , (1 - x)%) = 2(* x , ¥ t ) + (¥, F). (2.22) 



Next, 



2(* t ,* :ct ) + 2a||* t || 2 = ||* t || 2 r 1 +2a||v|/ t ||2 
-±d t (j|(l - z)* x || 2 + H-igi 2 + 26((1 - a:)*., (1 - as)*„)^ + (tf t) F). 



Next, 



(2.23) 



2(*«, + 2a(tf xi * t ) = fltH^H 2 + 2a(tf x , tf t ) = (^ x , ((1 - x) 2 V x ) x ) + V yy ) 

+bU x ,{{l-x) 2 * x ) y \ + b(* x ,((l-x) 2 * y ) x ) + {^ X ,F). (2.24) 



Now, 



9 X , ((1 - xf^ x ) x ) = -^,2(l-x)* x J + (1 - a;) 2 *.* 

= -(* x , 2(1 (((l-z) a **)*,*«) -|l*«ll a r 0J 
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((1 - xf* x ) x ) = - 1 (1 - x)* x ) - h* x \\ 2 r . 



Also, 



(**,tf w ) = ~(^x V .^y) = -T;\\^y\\ 2 r 1 , 



I 

ft^.ui-x) 2 *,)^ = -ftua-x) 2 ^,^ =o, 

((1 - ^) 2 ^)x) = -2&f* x , (1 - x)*y\ + bU x , (1 - x)H xy 
= -2b((l-x)V x ,$ ^ 



Therefore (|2.24l) becomes 

d t \\* x \\ 2 + (* x ,(l-x)*^) +i||^|| 2 ro + I||^ y || 2 ri +26^(1 -x)^ x ^ y ^j = -2a(«£ x ,V t ) + {-S> X: F). (2.25) 

All the boundary terms have the right sign to enhance the estimates. Therefore we ignore them. (See the Remark 
below.) Adding the simplified estimates ([2722]) . (|2T23]) . (j2~25j) gives 



d t ^a\\n 2 + ll^ll 2 + \Qj + Q+ \* x , (1 - x)^ a 
-2U[ (1 - x)* x , tf^ + 2(1 - o)(* B , * t ) - 2o||* t || 2 + (* + ** + *„, F) 

< constY||* x || 2 + ||*,|| 2 + ||vl/|| 2 + ||F || 2 V (2.26) 



where 

Q = ||(1- x)^|| 2 + ||*. s || 2 + 26^(1 - z)*x, (1 - •?)*< 
Since |6| < 1, there is a S > such that 

Q> 5(11(1 -z)*,|| 2 +||vM 2 ). 
Therefore (12.26)) gives us the required estimate for the energy norm 

E = a\m 2 + \\^4 2 + \q- 

Remark: If the boundary term in (|2.22[) had not been ignored then we would have obtained a stronger estimate for 
the energy 

F = F + ||*|| 2 ri 

which also controls the growth of \& on the boundary T\. 

We shall now prove that the problem is stable against lower order perturbations. We add an expression 

P = A^ x + B^ y + C^t + L>* 

to (j2~2Tj) . Then the estimates for (ET2"2]) , (j2~2Tt)) and (EHM)) will be changed by lower order terms 

(*, a^ x ) + (*, B^> v ) + (*, c* t ) + (*, mi) 

(*t, A* x ) + (* t , B^ y ) + (* t , (7* t ) + (* t , 
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Clearly, there is an energy estimate, provided we choose 2a > \C\. Thus the strip problem (|2.21|) is well posed. 
Now we start with 

2$ x4 = ((1 - x)H x ) x + $ yy + b ((1 - x)H x ) y + b((l~ x) 2 %) x + S(t, x, y). (2.27) 
We make the change of variables 

$ = e ax ^>, 

i.e., 

= e ax V x + ae ax V, <S> XX = e ax V xx + 2ae ax V x + a 2 e ax V, 
and set F = e~ ax S. Then we obtain (|2.21j) which is modified by 1Z 

2{^ x + a*) t = ((1 - xf^ x ) x + V yy + b (((1 - x)H x ) y + ((1 - xf^ y ) x ) +F + TZ. (2.28) 

Here 1Z consists of lower order terms, 

K = 2a(l - x) 2 ^ x + (a 2 (l - xf - 2a(l - x)) * + 2a6(l - xf^ y . 

Since (|2.21|) is stable against lower order terms there is an energy estimate for ()2.28|) . In the same way as in Sec. Ill Bl 
we can estimate the higher derivatives. This allows us to extend well-posedness to the variable coefficient problem 
and, locally in time, to the quasilinear problem. 

III. THE QUASILINEAR WAVE EQUATION ON AN ASYMPTOTICALLY FLAT BACKGROUND 

We now treat the null-timelike initial-boundary problem (|1.16[) for the quasilinear wave equation. We compactify 
the domain R < r < oo by the transformation x = 1 — R/r to obtain a strip problem < x < 1 with future null 
infinity T + at the boundary x = 1. In terms of the rescaled field $ = r$, the wave equation transforms into 

2d u d x i> - R- x d x (W(l - x) 2 d x $>) + (1 - x)R- x {d x W)$ 
+R- 2 D A ((1 - x) 2 W A d x $>) + R- 2 d x ((l - x) 2 W A D A $>) - R- 2 (l - x){D A W A )§ 

-R- 1 D A {e 2 f ) D A ^>) = -r 3 i?- x e^5. (3.1) 

Here we use the conformal 2-metric h A s and its inverse h AB to raise and lower indices of tensor fields on the spacelike 
(u = const., r — const.) spherical cross-sections . Up to lower order terms, (I3.1[) is a 3(spatial)-dimensional version 
of (|2.27l) where the y-coordinate has been replaced by the x^-coordinate on the spherical cross-sections and the 
i-coordinate has been replaced by the M-coordinate. In order for our treatment to apply to the quasilinear case, 
we assume that the metric coefficients (W, (3, W A , H A b) depend smoothly upon (4>, u, r, x A ) and that the source S 
depends smoothly upon ($, 9 a $, u, r, x A ), with non-singular Lorentzian geometry in the neighborhood of the initial 
data. 

We treat the modified problem resulting from the transformation $ = e ax ^. The same argument used in Sec. Ill Cl 
shows that this problem is stable with respect to lower order terms. We ignore these terms and thus obtain the strip 
problem 

2d u (d x V + a*) = R^d^Wil - x) 2 d x ^) - R- 2 D A ((1 - xfW A d x V) - R- 2 d x {{\ - x) 2 W A D A V) 

+R- 1 D A {e 2() D A ^) +F, %(0,x,x A ) = /, 0, x A ) = q , (3.2) 

where F = —r 3 R~ 1 e 2 Pe ax S. In order to treat (|3.2I) . we require that the physical space source has asymptotic behavior 
S = 0(r -3 ) so that F is square integrable over the strip. No boundary condition is allowed at the outer boundary 
Ti at x = 1 since I + is an ingoing characteristic surface. 

We obtain the required estimates for (|3. 21) by the same method used in Sec. Ill CI The data at the inner boundary 
T at x = may be absorbed into F so it suffices to treat the case q — 0. We define the inner product 

(*i,* 2 )= J dx j> dw*i* 2 
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and L2 norm ||^|| 2 = \l r ), where dw is the area element on the unit sphere. We write 

\\V A f = (V A ,V A ) = (h AB V A ,V B ). 

Since the spherical cross-sections are spacelike, their intrinsic 2-metric Iiab is positive definite so that || Va\\ serves as 
an Li norm for the angular components. We also need a metric norm for spacelike 3- vectors. In the standard Cauchy 
problem this is supplied by the intrinsic 3-metric of the spacelike Cauchy hypersurfaces. Since the characteristic 
hypersurfaces have a degenerate 3-metric, we take a different approach. We use the projection operator 7r£ = S^—t a dbU, 
where t a d a = d u , to define a 3-metric 7 a6 = n^n^g^. For the Bondi-Sachs metric (|1.15[) . the resulting components in 
the (it,r, x A ) coordinates are 
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-0, Y r = e- 2P W, Y A = ~e~ 20 r~ 2 W A , j AB = r - 2 h J 



Denoting x % — (r, X A ), this implies that 7 1 -' — > e 4 - 7 as r — > 00, where e 4 - 7 is the Euclidean 3-metric expressed in standard 
spherical coordinates. In the compactified coordinates, x l = (x, x A ) it is more useful to deal with the rescaled 3-metric 



~ ab 



1 



_ e 2/3 r 2^ab w ] 1 i c i 1 components 



7 a "=0, 7 XX = (1 - x) 2 W, ^ xA = -R- l (l - x) 2 W A , j AB = e 20 h AB . 

We then define 

\\V i t = (rf*V i ,V j ) 
which serves as an L2 norm for the x % components. Thus 

||d 4 *|| 2 = 11^/2(1 _ X )d x ^f + ||e^D A *|| 2 - 2R- 1 (V - x)d x V, (1 - x)W A D A 9 j . ( :-5.3) 

We also define the corresponding inner products and norms on the boundaries, e.g. 

(*i,*2)r = j> do;*i* 2 , ||*||r = (*,*) 



r- 



Because the radial coordinate r used in the Bondi-Sachs metric (|1.15p is a surface area coordinate, the conformally 
rescaled 2-metric h A B — t 2 9ab of the spherical cross-sections has determinant det(h A s) — det(q A B), where q A B is 
the unit sphere metric. Consequently, 

(V u D A D A y 2 ) = -(D A ^ U D A ^ 2 ) (3.4) 

and 

(V A ,D A *) = -(D A V A ,V) (3.5) 

where V (u, x, x A ) is any smooth vector field on the spherical cross-sections. These identities allow the necessary 
integration by parts. 

We derive the required estimates by freezing the dependence of the metric coefficients on (^,u, x) but we retain 
their dependence on x A so that W A and h A B remain smooth vector and tensor fields on the spherical cross-sections. 
We follow the procedure in Sec. Ill CI First, 

2(*,5 u ^*) + 2a(*,a u *) = -2(d x ^,d u 'f)+d u \\n 2 r 1 +ad u |tf|| 2 
= -R- 1 [W(l - x)d x ^, (1 - x)d x *) - R-^Da^W 2 + 2R- 2 ((1 - x)W A D A V, (1 - x)d x ^ ) + 



i.e. 

a„||*|| 2 ri +aa ll ||*|| 2 + i?- 1 ||a l *|| 2 = 2(^*,a u *) + (*,F). (3.6) 

Next, 

2(d u ^,d u d x ^) + 2a\\d u n 2 = \\dun 2 r 1 +2a\\d u n 2 
= -iR^du ( ||ir/ 2 (l - x)d x n 2 + \W 3 D A n 2 - 2i?- x ((l - x)d m V, (1 - x)W A D A A + {d u V,F) , 
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so that 



Next, 



|3 u *|| 2 ri +2 a ||0 u *|| 2 + ^R-'duWd^W 2 = (d u *, F). (3.7) 



2(3**, d u d x ^) + 2a(d x V, 8*9) = d u \\d x V\\ 2 + 2a(d x 9, d u 9) 
= R- 1 (d x ¥, d x (W(l - x) 2 d x 9)\ + iT 1 (d x 9, D A (e 2fj D A 9) 

-R- 2 (tf x , D A ((1 - x) 2 W A d x 9)) - R- 2 (d x 9, d x ((l - x) 2 W A D A 9) ) + (fit*, F). (3.8) 



As shown in in Sec. Ill CI 



d x 9, d x (W(l - x) 2 d x 9)\ = - (d x 9, W(l - x)d x *) - ^\\W^ 2 d x n 2 r - 



Also, 



(d x 9, D A (e 2p D A 9)) = ~\\e^D A 9\\ 2 ri (3.9) 

and 

/ d x *,D A ((l-x) 2 W A d x V) S ) + (d x 9,d x ((l-x) 2 W A D A 9)\ = -2^(1 -x)0 x 9,W A D A W 

Therefore Q3.8P becomes 

fi»R*|| 2 + R- 1 (d x 9, W(l - x)d x *\ + iR-'WW^d^fro + iR-^Dj&fr, 

-2R~ 2 ((l-x)d x 9,W A D A 9\ = -2a(d x 9,d u 9) + (d x 9,F). (3.10) 

As before, the boundary terms have the right sign to enhance the estimates so that we can ignore them. Adding 
the simplified estimates (13.61) . p.7[) , (|3.10p gives 

d u (a\\n 2 + \\d x n 2 + jiT 1 1| 9;* + R^d^f + R' 1 (dx*, W(l - x)d x V 
= 2R- 2 ({I - x)d x 9, W A D A *) + 2(1 - a)(d m % d u 9) - 2a||d u *|| 2 + (* + U 9 + d x 9, F) 



< const. I ||<9 X *|| 2 + \\D A 9\\ 2 + ||*|| 2 + \\F\\ 2 j . (3.11) 

Therefore (I3.11[) gives us an energy estimate provided that the 3-metric 7 y has (+ + +) signature, so that ||9i*|| is 
a norm for the gradient fij* = (d r 9,d A 9). This is equivalent to the requirement that the principal part of the wave 
operator reduce to an elliptic operator in the stationary case where the it-derivatives vanish. Since 7 IJ is asymptotic 
to the Euclidean metric as r — > 00, this positive-definite condition is satisfied throughout some exterior domain. 

Estimates for the higher derivatives of 9 and stability against lower order perturbations follow from the same 
arguments given in Sec. [TTJ This establishes the well-posedness of the worldtube-nullcone problem for the case of 
smooth variable coefficients. The extension of well-posedness, locally in time, for the quasilinear case then follows 
from the standard techniques referred to in Sec. [TT] . 

For a mass M Schwarzschild geometry, 7 rr = e~ 2 ^W = 1 — 2M/r so that positive-definiteness of the 3-metric ^ 
breaks down at r = 2M where the worldtube becomes null. In this limiting case of the double-null problem, the 
<9 U 1 1 3^ "J 1 1 2 term in (|3.11[) suffices to provide the required estimate. However, for R < 2M the "worldtube" is spacelike 
and the problem must be treated differently. 
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